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We report on measurements and modeling of the mode structure of tunable Fabry-Perot optical 
microcavities with imperfect mirrors. We find that non-spherical mirror shape and finite mirror 
size lead to loss, mode deformation, and shifted resonance frequencies at particular mirror separa¬ 
tions. For small mirror diameters, the useful cavity length is limited to values significantly below 
the expected stability range. We explain the observations by resonant coupling between different 
transverse modes of the cavity and mode-dependent diffraction loss. A model based on resonant 
state expansion that takes into account the measured mirror profile can reproduce the measurements 
and identify the parameter regime where detrimental effects of mode mixing are avoided. 

PACS numbers: 42.25.Fx, 42.15.Eq, 42.50.Wk, 42.60.Da 


Fabry-Perot optical microcavities built from micro- 
machined concave mirrors offer a powerful com¬ 

bination of small mode cross section, high finesse, and 
open access. This has proven to be beneficial for experi¬ 
ments covering a broad range of topics, including cavity 
quantum electrodynamics with cold atoms mm , ions 
nnnu, and solid-state-based emitters [IlHiH], as well 
as cavity optomechanics [TOHSI] and scanning cavity mi¬ 
croscopy [22]. Various techniques have been developed 
to produce concave, near-spherical profiles as mirror sub¬ 
strates, including CO 2 laser machining El El [a na ng, 
chemical etching [31 El] , focused ion beam milling laES], 
and thermal reflow EHEI]- A small cavity mode cross 
section is achieved by realizing microscopic surface pro¬ 
files with radii of curvature Tc ^ 5 — 500 pm and pro¬ 
file diameters typically a factor 2 — 10 smaller. In this 
regime, the extension of the cavity mode can be compara¬ 
ble to the effective mirror diameter, and the finite mirror 
size becomes relevant. Furthermore, the different fabri¬ 
cation processes typically yield profiles that deviate from 
a spherical shape, and excessive surface roughness may 
be present. In addition, coating defects and particles on 
the mirror surface can disturb cavity performance under 
real conditions. Overall, the mode structure of open- 
access microcavities will be affected by the details of the 
mirrors, and an in-depth understanding of the relation 
between mirror imperfections and cavity performance is 
required for the successful application and the potential 
improvement of such resonators. 

In this work we study the consequences of finite 
mirror size and non-ideal shape on the performance 
of laser-machined, fiber-based Fabry-Perot microcavities 
[Tj. Their mirrors are characterized by surface profiles 
with low microroughness in the range of 1 — 2 A, a near- 
spherical central part, and an overall shape that is well 
approximated by a Gaussian. 

We perform measurements of the cavity transmission 


and finesse across the entire stability range for several 
cavities. We find that for short mirror separation, the 
cavities are mostly immune to mirror imperfections, and 
the fundamental cavity mode closely resembles a Gaus¬ 
sian mode. However, at particular mirror separations, 
the cavity shows a significantly reduced finesse, and the 
performance depends for example on the precise laser 
wavelength. Furthermore, for small mirror size, we ob¬ 
serve that the distance range where the finesse remains 
high is significantly smaller than the stability range ex¬ 
pected from the mirror radius of curvature. 

We accurately reproduce these observations with a 
model based on resonant state expansion [28ll3Q] . where 
we take into account the measured mirror profile. The 
model shows that the observed behavior can be consis¬ 
tently explained by (near-) resonant coupling between 
different transverse modes of the cavity, caused by the 
non-ideal shape and finite size of the mirrors. The ad¬ 
mixture of higher order modes, which suffer from diffrac¬ 
tion loss due to their larger size, introduces loss to the 
fundamental mode. Based on these results, we identify 
the parameter regime for Gaussian-shaped mirrors where 
detrimental effects of mode mixing remain negligible. 


FIBER-BASED MICROCAVITY 

The cavity design is depicted schematically in Fig. 
(a): The resonator consists of a curved micromirror ma¬ 
chined on the end-facet of a single mode optical fiber and 
a macroscopic plane mirror. Both the commercial plane 
mirror substrate and the fiber surface are coated with a 
highly reflective dielectric coating for a center wavelength 
of 780 nm, where a finesse of T ^ 60000 can be reached. 
In this configuration, the planar mirror serves as a near¬ 
ideal reference mirror, which permits us to study the ef¬ 
fects of the micromirror alone. The light of a grating- 
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FIG. 1: (a) Schematic setup showing the machined fiber, the plane mirror, and the optical mode, (b) Mirror profile as 

measured with a white-light interferometer. Dashed circle: typical mode size (1/e^-diameter) on the mirror, (c) Profile sections 
(as indicated with dashed gray lines in (b)). Red (Orange): cut along x (y) -direction. Dashed lines: Gaussian fit. (d) Residuals 
of the fits in (c). (e) Residual of a 2D parabolic fit to the region of the profile typically covered by the cavity mode. 


stabilized diode laser is coupled into the cavity through 
the fiber, and light transmitted through the plane mir¬ 
ror is collimated and detected with an avalanche photo 
diode. The whole stability range of the cavity can be 
covered with sub-nanometer resolution with a piezo step 
drive linear positioning stage (PI LPS-24), onto which 
the fiber is mounted. In addition, a shear piezo crystal 
is used for scanning the cavity length over the resonance. 
The plane mirror can be laterally scanned with an XY 
piezo table (PI P-541.2SL) over one hundred microme¬ 
ters. A mirror mount allows for angular alignment of the 
cavity. 

A white-light interferometric image of the laser- 
machined depression on the fiber surface is shown in 
Fig. 0 (b). The image is taken with a home-built in¬ 
strument with a lateral resolution of 560 nm and a ver¬ 
tical resolution of 0.1 nm (rms). Cuts through the pro¬ 
file along the dashed gray lines are shown in Fig. 0(c). 
The dashed white circle illustrates the 1/e^-diameter of 
the fundamental mode for a mirror separation d ~ rc/4. 
A Gaussian fit to the surface data (lines in Fig. 0(c)) 
shows good agreement, and only nanometer-scale devia¬ 
tions can be seen from the residuals (Fig. (d)). The 
central part of the profile can be well approximated by 
a parabola. Figure (e) shows the residual of a two- 
dimensional parabolic fit to the data. Certain localized 
imperfections are present on the fiber surface, as well as 
an overall shape deviation. 

Figure [^(c) also illustrates that the profile is not rota- 
tionally symmetric but rather has elliptical contour lines. 
This leads to Hermite-Gauss modes to closely resemble 
the eigenmodes of the cavity, and to a splitting of higher- 
order transverse modes of equal order. Additionally, the 
ellipticity splits each cavity resonance into a linear po¬ 
larization doublet [23]. We use polarization optics be¬ 
fore the fiber to select one of the modes for evaluation. 
The surface shown in Fig. (b) is rotated such that the 


principle axes of the profile coincide with the coordinate 
axes. The minimal radius of curvature in the center is 
found to be = 161 pm in x— and = 201 pm in 
^—direction. 

EXPERIMENTAL RESULTS 

We study the cavity performance by measuring the fi¬ 
nesse for each accessible axial mode order. To ensure that 
local variations of the mirror coating of the planar mirror 
do not influence the result, we determine the finesse at 
25 different positions within an area of 30 x 30 pm on 
the planar mirror and evaluate the most prevalent value. 
Diffraction loss, which arises as result of mode mixing, 
leads to a decreased finesse according to 


where T denotes the total transmission of the mirrors, A 
the total absorption loss, and D the diffraction loss due 
to the micromirror. In our experiment T +A ^ 100 ppm. 

A typical measurement of the finesse of the fundamen¬ 
tal mode as a function of mirror separation is shown 
in Fig. 2. In such measurements, we typically observe 
three different regimes: For small mirror separation d ^ 
rc/4, (axial mode number q < 130 for the measurement 
shown), the finesse stays approximately constant with 
only a slight overall decline. For intermediate mirror sep¬ 
arations, individual axial mode numbers q show large ad¬ 
ditional loss. For mirror separations d ^ rc/2, an abrupt 
drop of the cavity transmission and finesse is observed, 
with few moderately-working mode orders appearing for 
larger d. 

We observe that increased loss appears mainly when 
higher-order transverse modes become degenerate with 
the fundamental mode. Therefore it is instructive to 
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FIG. 2: Cavity finesse as a function of axial mode order de¬ 
termined from measurements at 25 different positions within 
an area of 30 x 30 pm on a large planar mirror and taking 
the most prevalent values. 


study cavity transmission spectra for different mirror sep¬ 
arations. A typical spectrum of a cavity at high finesse 
is shown in Fig. 3 (a), where we probe the cavity with a 
narrow-band laser and tune the mirror separation across 
one free spectral range. To map out the mode profile, we 
raster-scan a nanoparticle placed on the large mirror us¬ 
ing the cavity mode and evaluate the introduced loss from 
the cavity transmission of each of the modes [22]. The 
mode shapes clearly resemble those of Her mite-Gaussian 
modes orientated along the principle axes of the mirror 
profile. For a given mode with transverse mode order 
(m, n) and axial mode order q, the cavity resonance fre¬ 
quencies for a spherical mirror cavity are given by 
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where = arccos y 1 — d/r^'^^ is the Gouy phase. 

Note that the degeneracy of modes with the same trans¬ 
verse mode order m + n is lifted by the ellipticity of the 
profile, leading to families of modes with m + n +1 mem¬ 
bers. Figure 3 (b) shows spectra like the one in (a) as 
a function of cavity length for 856 nm probe light, where 
T ~ 1200 to improve the visibility of the resonances. The 
logarithmic color scale is set for each cavity length to op¬ 
timize the signal-to-background ratio. With increasing 
mirror separation, the spacings between the transverse 
modes increase and eventually, higher-order modes can 
become resonant with the next fundamental mode. This 
is given when the differential Gouy phase fulfills 


= jn 


(3) 


where j is an integer. 

When evaluating the measured transverse modes, we 
find that their frequencies deviate from the spectrum 


FIG. 3: (a) Gavity transmission spectrum as a function of rel¬ 
ative cavity detuning covering one free spectral range, where 
T 60.000. The insets show the mode functions as measured 
by scanning cavity microscopy, (b) Gavity transmission spec¬ 
tra as a function of axial mode order q for T ~ 1200. Loga¬ 
rithmic color scale. White: model for modes (0, m + n) and 
(m + n, 0). 


given by Eq. 2, and Eq. 3 fails to predict the positions 
of the observed resonances. An accurate description is 
possible by the model discussed in the next section. The 
white lines in Eig.[^(b) show the predicted resonances of 
the modes (0, m + n) and (m + n, 0) for the lowest few 
mode orders mn. 

A closer look at the crossover of modes 04 and 13 with 
the fundamental mode is taken as an example for the 
typical avoided crossing behavior found when mode cou¬ 
pling occurs. Eig. 4 (a) shows spectra covering the re¬ 
gion around the fundamental mode for every axial mode 
order from q = 198 to q = 220 for T ~ 1200. Close 
to resonance, the transmission of the fundamental mode 
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FIG. 4: (a) Measured cavity transmission spectra as a func¬ 

tion of relative cavity detuning. Avoided crossings between 
the fundamental mode and the fourth order modes 04 and 13 
are visible, (b) Maximal transmission of each of the spectra 
shown in (a), normalized to the largest value. 


decreases (see Fig. 4 (b)), while the transmission of the 
higher order mode increases until both have approxi¬ 
mately equal height and minimal separation at the point 
of resonant coupling. The coupling is accompanied with 
an increased linewidth and thus a reduced finesse. At the 
anticrossing between the 00 and 04 mode, we observe a 
mode splitting of 8.6 GHz. We have modeled the cou¬ 
pling of the particular mode pair (see below) and find 
a value of 8.6 GHz, matching the measurement within 
errors. 


MODELING 


matrix M, which accounts for the change a mode under¬ 
goes during one round trip through the cavity, the task 
reduces to an eigenvalue problem 


7 .^. = M^i. 


(5) 


The mixing matrix M^^k' = exp(— 4 i 7 r(i/A)l/c,/c/H/c^/c/ has 
elements given by mode overlap integrals taken over the 
finite extent of the micromirror. 




k,k' 


/ xo pyo 

-Xo J-yo 


(6) 


Here, (xo,^o) denote the extension of the mirror, ± in¬ 
dicates the sign of the phase factor of 4>/c, A{x,y) is the 
deviation of the mirror profile from a planar surface, and 
Zm is the location of the micromirror on the optical axis. 
We assume that the respective expression for the planar 
mirror is an identity matrix. 

Using the Hermite-Gauss modes for the expansion im¬ 
plies the paraxial approximation, where the isophase sur¬ 
face is parabolic (with some deviation due to the Gouy 
phase) rather than spherical. However, the paraxial ap¬ 
proximation does not hold at large separation from the 
optical axis where the two shapes differ. In fact, when in¬ 
cluding non-paraxial terms, one finds that a spherical ge¬ 
ometry is indeed the most desirable m- For a spherical 
mirror which covers an entire half-space, Bk^k' = ff/c,/c '5 
and the eigenvalues 7 ^ corresponding to the eigenmodes 
are unity. As soon as A{x^y) deviates from spheri¬ 
cal or xo^yo is finite, M has off-diagonal elements and 
transverse-mode mixing occurs. 

For an accurate treatment of our experiment, we use 
the measured surface profile (Fig. 1(b)) for A{x^y). To 
find a suitable basis set ^k for each mirror separation, we 
numerically maximize |Mo,o| by varying the mode waist 
wq of 4>oo. For a given mirror separation, the obtained 
optimal Wq corresponds to an effective radius of curvature 
^c,eff of ^ spherical mirror. The result for the profile 
investigated here is displayed in Fig. |^(c), showing that 
^c^eff larger than = 161 pm and increases with d. 
Consequently, the stability range is expected to extend 
beyond the limit of d = Tc- 

The diffraction loss Di of a mode can be directly 
obtained from the corresponding eigenvalue 7 ^, 


For cavity mirrors where the surface profile can be 
treated as a perturbation of a spherical shape, an effec¬ 
tive approach is to describe the real eigenmodes as a 
series expansion of Hermite-Gauss modes <!>/., 

^i = 'Y^Cik^k, k = {m,n). (4) 

k 

Following the approach of Kleckner et al. [29], we deter¬ 
mine the new eigenmodes and the corresponding reso¬ 
nance frequency and loss. Introducing the mode-mixing 


A = 1-|74^. (7) 

Inserting this into Eq. gives the finesse of this mode. 

The obtained finesse of the fundamental mode using 
the profile shown in Fig. 1 (b) is plotted in red in Fig. 
5 (a), where we use the measured mirror transmission 
and loss. For direct comparison, we show a measurement 
of the cavity finesse as obtained from the resonant cav¬ 
ity transmission (Tc oc J^^) at a wavelength of 780 nm. 
For each data point, we measure 22500 positions on a 
30 X 30 pm area of the plane mirror and take the most 
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FIG. 5: (a) Data points: Cavity finesse extracted from 22500 measurements for each datapoint. Red: Simulation of the cavity 

finesse as a function of the mirror separation for a measured surface profile assuming mirror transmission and absorption to be 
100 ppm. Green: Simulation for a Gaussian prohle fitted to the measured surface, (b) Composition of the ground mode from 
Hermite-Gaussian modes contributions |cofc|^ are shown. Black: |coo|^- Black numbers denote the transverse mode order 
of other important contributions. Shaded areas show impact zones of specihed mode orders. Inset: |^o|^ for a region of low 
coupling {q = 60). (c) Effective radius of curvature sensed by the fundamental cavity mode as a function of cavity length, 
(d) Examples of mode shapes of the fundamental mode |^o|^ for selected mirror separations indicated by Roman numbers in 
(b) exhibiting large coupling. The edge length is 40 pm . 


prevalent value. The rise of the finesse for short cavities 
can be attributed to a systematic error of the measure¬ 
ment. The overall shape, the exact position of localized 
finesse dips, and the decrease around q = 230 can be 
reproduced by the simulation with a high level of de¬ 
tail. However, to match the data, the lateral size of the 
mirror profile had to be rescaled by about 2.5% for the 
simulation. The same correction has been made for the 
simulation of the spectrum shown in Fig. 3, where the 
normalized phase of the eigenvalues 7 ^ is plotted. The 
mismatch might result from a calibration uncertainty of 
the interferometric surface reconstruction. In contrast, 
non-paraxial corrections of the mode frequencies [32] are 
estimated to lead to relative frequency changes below 
10 “^ for the current cavity geometry. 

The computed eigenvectors contain information about 
the composition of the system’s eigenmodes from 
Hermite-Gaussian modes according to Eq. In Fig. 5 
(b), the coefficients |co/cP giving the contributions to the 
fundamental mode are plotted as a function of mirror sep¬ 
aration. The Gaussian mode $00 is clearly the dominant 
one, and for most cavity lengths, the ground mode shows 
negligible deviation from it (see inset in Fig. |^(b)). How¬ 
ever, for certain distinct mirror separations where reso¬ 
nant mode mixing occurs, higher order modes can have 


significant contributions and lead to a severe distortion 
of the fundamental mode (Fig. 5(d)). The larger spatial 
extent of higher-order modes with Wk ~ woVk + 1 causes 
larger diffraction loss, from which also the fundamental 
mode suffers under coupling conditions. Notably, the lo¬ 
cations of high loss and strong mode mixing do mostly 
but not necessarily coincide (see below and [29]). 

Regions of impact of certain mode families (shaded ar¬ 
eas) cover a significant fraction of the stability range. 
Still, for applications where the exact mirror separation 
is not essential, extended regions of negligible mode mix¬ 
ing remain. The different influence of e.g. mode orders 
4 and 5 can be attributed to the larger values of Ho,/c 
for even modes due to the symmetric mirror profile and 
modes with smaller mode index. Also, for larger mode in¬ 
dex differences, the differential Gouy phase evolves faster 
and the resonance condition is sharper. The coupling 
strength can be directly inferred from the mode splitting 
at an avoided crossing. 

It is instructive to compare the results with a calcu¬ 
lation for a profile obtained from a Gaussian fit to the 
measured fiber surface (green solid line in Figj^ (a)). 
The smooth surface does not lead to the overall weak 
decline for increasing mirror separation, and the sharp 
features at intermediate d are missing. Yet, the finesse 
























































6 



FIG. 6: Simulated finesse for different mirror profile size 
at fixed radius of curvature Vc using various ratios of Wc/a 
defined at d = rc/2. Solid lines: Gaussian profile. Dashed 
lines: parabolic profile. We assume additional loss of tt x 10“^ 
for each mirror. 

decrease around d = rc/2, which effectively limits the 
stability range, is reproduced. The difference can be ex¬ 
plained by the presence of additional (and in particular 
asymmetric) surface deviations with mostly larger spa¬ 
tial frequencies and particle-like features in the measured 
profile. High spatial frequencies couple the fundamental 
mode to many transverse modes with large mode index, 
causing a smooth finesse decrease and significant reso¬ 
nant mixing for particular modes. 

The observed behavior is not limited to the particu¬ 
lar parameters used in our experiment, but is a general 
property related to the profile shape and size. Consid¬ 
ering a Gaussian profile with 1/e radius a, depth t, and 
Vc = l(2t\ and assuming a cavity with d = rcl‘^ where 

the mode radius on the curved mirror is Wc = 
one finds that the relative mode size depends only on the 
profile depth for a given wavelength, Wcjd = 

We perform simulations for a profile with fixed Tc and 

ellipticity e = y 1 — jr^c^ = 0.26 [23] and vary rCc/a. 
The resulting finesse is shown in Fig. [^ While profiles 
as small as a = 2wc already achieve performance not 
limited by diffraction for small mirror separation, it re¬ 
quires a profile radius a > Awc to extend this range to 
d = rd^ and a > 10 Wc to avoid mode mixing over the 
entire stability range. For comparison, we also perform 
simulations for a rotationally symmetric parabolic profile 
with Vc = Vc ^ and an edge length of 2a. While the overall 
behavior is similar, the calculation for Wc = 0.56 a shows 
that resonant transverse-mode mixing can also lead to a 
reduction of diffraction loss [29| . This can be understood 
by the destructive interference between the fundamental 
and the higher-order mode at the outer part of the mode, 
reducing the effective mode size [22]. 

For obtaining the presented data, modes up to order 
m + n = 20 are included in the calculation. Using more 


modes does not significantly alter the results and strongly 
increases the computation time, which grows approxi¬ 
mately as (m + n)^. For a fiber profile of 400 x 400 
pixels, the simulation of 300 mirror separations could be 
conducted with a personal computer within a few hours 

m- 


CONCLUSION 

Our results comprise an extensive analysis of the per¬ 
formance of optical cavities with non-ideal mirrors. We 
have shown that localized reduction in cavity finesse, fre¬ 
quency shifts, and mode shape distortions are the con¬ 
sequences of non-ideal mirror shape and finite size. The 
behavior can be accurately modeled with a mode expan¬ 
sion method when using the measured mirror profile as 
an input. The demonstrated approach provides a pow¬ 
erful tool for analyzing a given cavity geometry and for 
predicting cavity performance without the need for de¬ 
tailed experimental characterization. This is particularly 
important for experiments where minimal mode volume 
and ultimately small mirror profiles are desired, as well as 
for cold atom, ion trap, and cavity optomechanics exper¬ 
iments, where larger mirror separations in combination 
with small radii of curvature are beneficial. The calcu¬ 
lations provide improved accuracy for the determination 
of emitter-cavity coupling strength as well as detailed in¬ 
formation about possible sample-induced scattering and 
loss. Finally, the approach offers an efficient route for the 
design of novel cavity geometries with non-trivial prop¬ 
erties, such as single-transverse-mode operation [33] or 
mode imaging [34]. This opens the potential for light 
modes to be individually tailored for specific applications. 
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